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ABSTRACT 

We examine the infrared properties of 43 high rcdshift (0.1<z<1.2), infrared- 
luminous galaxies in the Extended Groth Strip (EGS), selected by a deep 70 ^m 
survey with the Multiband Imaging Photometer on Spitzer (MIPS). In addition and 
with reference to starburst-type Spectral Energy Distributions (SEDs), we derive a 
set of equations for estimating the total infrared luminosity {Lm) in the range 8- 
1000 /zm using photometry from at least one MIPS band. 42 out of 43 of our sources' 
optical/infrared SEDs {Xobserved < 160 /im) are starburst-type, with only one object 
displaying a prominent power-law near-infrared continuum. For a quantitative analy- 
sis, models of radiation transfer in dusty media are fit onto the infrared photometry, 
revealing that the majority of galaxies are represented by high extinction, A^, > 35 
and for a large fraction (~ 50 per cent) the SED turns over into the Rayleigh- Jeans 
regime at wavelengths longward of 90 /im. For comparison, we also fit semi-empirical 
templates based on local galaxy data, however, these underestimate the far-infrared 
SED shape by a factor of at least 2 and in extreme cases up to 10 for the majority 
(~ 70 per cent) of the sources. Further investigation of SED characteristics reveals that 
the mid-infrared (70/24 /xm) continuum slope is decoupled from various galaxy prop- 
erties such as the total infrared luminosity and far-infrared peak, quantified by the 
L160/L70 ratio. In view of these results, we propose that these high-redshift galaxies 
have different properties to their local counterparts, in the sense that large amounts 
of dust cause heavy obscuration and are responsible for an additional cold emissive 
component, appearing as a far-infrared excess in their SEDs. 

Key words: galaxies: general galaxies: high-redshift galaxies: starburst galaxies: 
photometry 



1 INTRODUCTION 

The classes of galaxies known as Luminous, Ultraluminous 
and Hyperluminous Infrared Galaxies (LIRGs, ULIRGs and 
HyLIRGs), have been a major subject of focus in observa- 
tional cosmology since initial, groundbreaking studies such 
as those of Soifer et al. (1986) and Sanders et al. (1987). 
These objects, discovered in abundance by the Infrared As- 
tronomical Satellite (IRAS), were found to possess an array 
of extreme properties, including an excessive energy out- 
put, effectively manifested as a large infrared/optical ratio 
(Sanders et al. 1989) . It is now widely accepted that the com- 
bination of a powerful UV-photon source and the presence 
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of substantial amounts of dust are responsible for such ex- 
treme luminosities of the order of 10^^-10^^ Z/0 (e.g. Sanders 
& Mirabel 1996, Genzel et al. 1998). Determining the na- 
ture of the central energy source has been the subject of 
numerous subsequent studies, very often supplemented by 
multiwavelength data (e.g. Gregorich et al. 1995, Genzel et 
al. 1998, Rigopoulou et al. 1999, Klaas et al. 2001, Tacconi 
et al. 2002, Alonso-Herrero et al. 2005). These have revealed 
populations dominated by powerful starburst activity, but 
with a non-negligible Active Galactic Nucleus (AGN) con- 
tribution found to increase with bolometric luminosity (e.g. 
Veilleux et al. 1997, Tran et al. 2001, Brand et al. 2006). 

The processes responsible for energy production and 
radiation transfer in the interstellar medium of such ob- 
scured systems are primarily evaluated by studying galaxy 
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Spectral Energy Distributions (SEDs), sometimes replac- 
ing spectroscopic diagnostics in determining the nature of 
the central engine. For infrared-luminous objects, these pro- 
cesses are greatly influenced by the composition, tempera- 
ture and distribution of interstellar dust. Quantifying the 
infrared energy budget is, therefore, key in determining the 
generic properties of a galaxy and, as it is directly cou- 
pled to the central energy source, important in estimating 
the integrated cosmic star formation history. Cases where 
photometric data is scarce, benefit from various methods 
of data extrapolation and interpolation: construction of in- 
frared SEDs has been a popular approach of modelling emis- 
sion from dusty systems, either directly from first principles 
or semi-empirically, bridging together contributions from 
different parts of the interstellar medium (e.g. Guiderdoni et 
al. 1998, Calzetti et al. 2000, Rowan- Robinson 2000). Such 
models have enabled the examination of sources with ex- 
treme properties and pronounced infrared luminosity, either 
by solely considering star formation processes or including 
contribution from AGN. In addition, establishing correla- 
tions between monochromatic flux densities or luminosities 
has proven extremely advantageous and has conveniently 
been used as a first order approximation of the infrared en- 
ergy budget, as well as a diagnostic for the physical pro- 
cesses in galaxies, adopted extensively in evolution studies 
(e.g. Franceschini et al. 2001). 

In this paper we investigate the properties of 43 
infrared-luminous galaxies from one of the first deep Spitzer 
far-IR surveys. This sub-sample is part of a population of 
178 sources detected at 70 /im by the Multiband Imaging 
Photometer on Spitzer (MIPS) (hereafter the 70 ^m pop- 
ulation), down to a limiting flux density of 4mJy at 5<j. 
In Symeonidis et al. 2007 (hereafter S07), we gave a brief 
overview of the sources' properties and derived Star Forma- 
tion Rate (SFR) estimates. Here, we aim to gain qualitative 
insight into their physical properties, primarily with respect 
to dust, comparing them to similar objects in the local uni- 
verse. For this purpose and to quantitatively characterise 
the sample, we fit our photometry with three starburst dust 
models and a set of empirical local galaxy templates. 

The paper is set out as follows: Section 2 introduces our 
selection criteria for the sub-sample, including an overview 
of infrared colours. Section 3 analyses and compares the 
model SED templates, associating the quality of fits to the 
properties of the sources. We investigate the infrared energy 
budget, with reference to SED characteristics such as the 
continuum slope and far-IR peak, in section 4. Finally, sec- 
tion 5 focuses on the derivation of a set of equations to cal- 
culate the total infrared luminosity in the range 8-1000 /im, 
with at least one MIPS band, for starburst-type sources. 
Our conclusions and discussion are presented in section 6. 
In all subsequent calculations, we have employed the follow- 
ing values: Ho = 71 kms~^Mpc~^, Qm = 0.3 and JIa ~ 0.7. 



2 INTRODUCING THE SAMPLE 

2.1 Data sets and Selection Criteria 

Our main data set is a 70 nm image of the Extended Groth 
Strip (EGS) field, taken by MIPS. In this -0.5 deg^ field, 
178 sources were recovered down to a limiting flux density 



of 4 mjy at 5a. We direct the reader to Davis et al. (2007) 
and S07 for a detailed description of all data sets acquired 
through the All- Wavelength Extended Groth Strip Interna- 
tional Survey (AEGIS). In this section, we briefly refer to the 
Deep Extragalactic Evolutionary Probe 2 (DEEP2) spectro- 
scopic survey, the Spitzer Infrared Array Camera (IRAC) 
survey and the MIPS survey. 

70 and 90 per cent of the MIPS 70 /im area is also cov- 
ered by the IRAC survey and MIPS 24 and 160 fim images, 
respectively. Within the overlap areas of all 8, 24 and 160 /xm 
images, we recover 100 per cent of our targets, a total of 122 
sources (hereafter, the full sample). Of these, 103 also appear 
in the remaining three IRAC bands (3.6, 4.5 and 5.8 /im). 
DEEP2 obtained optical spectra for ^^14,000 objects in the 
EGS with a limiting magnitude of Rab = 24.1, using the 
DEIMOS spectrograph on the Keck telescopes (Davis et al. 
2003). In the overlapping area (~75 per cent) between the 
MIPS 70 Atm image and the DEEP2 EGS survey, 43 rehable 
(quality factor ^3) spectroscopic redshifts (zspsc) were re- 
trieved in the range 0.1 < z < 1.2, with a mean of 0.57 and 
a median of 0.51. For the remaining objects, some made the 
photometric cut but were assigned a poor quality redshift. 

For the purpose of this work we select a sub-sample 
of the 70 fim population (table 1) (hereafter, the redshift 
sample), which we require to be within the common area 
of the IRAC 8 fim and the MIPS 24 and 160 fim images, 
in order to sample the entire available infrared SED. The 
criteria are as follows: 

• 8 /im photometry from the IRAC survey 

• 24 and 160 /xm photometry from the MIPS surveys 

• a reliable spectroscopic redshift from DEEP2 

It is the latter criterion which forms the basis of the se- 
lection, as the former two are always satisfled within the 
overlap regions. 34 objects in the redshift sample also ap- 
pear in the 3.6, 4.5 and 5.8 fim IRAC bands. The MIPS (24, 
70 and 160 fim) colours of the redshift and full samples are 
similar (figure [1]); apart from a few objects at the extremes 
of the 160/70 colour, for the majority of the population the 
received flux at 160 fim is more than twice that at 70 
and shows no dependence on redshift or the strength of the 
24 /im mid-IR emission. 

2.2 Spectral Energy Distributions 

The SEDs of the 43 sources in the redshift sample are shown 
in figure [2] They are separated into three groups with re- 
spect to the total infrared luminosity in the 8-1000 range 
— Starbursts (10^" - IO^Lq), LIRGs (10^-10^^ Lq) and 
ULIRGs (10^^ - 10"Lq) (see section 4). The SED of an ob- 
ject which displays a clear power-law near-IR continuum is 
in black. This was eliminated through examination of IRAC 
colours commonly used to identify objects whose SED is 
AGN-dominated (e.g. Lacy et al. 2004, Stern et al. 2005). 

An obvious trend in the SEDs is the reduction of the 
stellar bump with respect to the infrared continuum, for ob- 
jects of increasing luminosity and, as this is a flux limited 
survey, redshift. A similar observation was made by Alonso- 
Herrero et al. (2006), who found the stellar bump greatly re- 
duced or completely wiped out in local ULIRGs. This could 
be a direct consequence of an increasing AGN versus star- 
burst dominance, indicating the tendency towards a power- 
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Table 1. The redshift sample (43 sources) - Object cardinal numbers with the prefix EGS70 
referring to the 70/im EGS survey, coordinates, spectroscopic redshifts, R-band magnitudes ( AB) , 
measured flux densities at 24, 70 and 160 /um. The MIPS 24 /xm calibration uncertainty is 10 per 
cent, whereas for MIPS 70 and 160 /im it is 20 and 30 per cent, respectively. 
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law near-IR continuum or simply due to the fact that the 
more infrarcd-lumiuous galaxies will be fainter in the opti- 
cal because of increased extinction. For the lower-redshift 
objects, the 6.2 and/or 7.7 jim PAH features are redshifted 
into the 8/im band, which displays a clear increase in flux. 
Contrary to the pronounced differences in the ncar-IR, the 
mid/far-IR (24-160 /im) part of the SEDs makes them indis- 
tinguishable; this has also been noted by other authors, e.g. 
Klaas et al. (2001) and it is potentially indicative of the fact 
that far-lR emission is independent from the nature of the 
central energy source. 



3 SED FITTING 

In past years, there have been numerous efforts to char- 
acterise and reproduce galaxy SEDs with model templates 
(e.g. Rowan-Robinson 1980, 1995, Appleton et al. 2004). 
These have been developed either semi-empirically or by ap- 
plying detailed radiation transfer calculations; in either case 
there has been reasonable success with respect to the mod- 
els' representation of real systems (e.g. Efstathiou Rowan- 
Robinson 1990, Granato 1994, Silva 1998, Rowan-Robinson 
et al. 2005, Efstathiou, Rowan-Robinson & Siebenmorgen 
2003, Siebenmorgen k Kriigcl 2007). Our motive for SED 
fitting is to determine the infrared energy budget of each 
source and gain insight into its physical properties. Since the 
redshift sample mostly displays SEDs typical of starburst 
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Figure 1. Infrared colours for the full sample - ratio of 160 /xm 
vfv to 70 /xm i/fi, versus 24 /im i^/jy. Asterisks represent R> 24.1 
sources with no Zspec , while the rest of the symbols are divided 
into z^O.5 (empty circles), z > 0.5 (filled circles) and R<24.1 
sources with no Zspec (crosses). The blue diamond symbol corre- 
sponds to the object we identified as AGN-dominated (see dis- 
cussion in section 2.2). 
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Figure 2. The optical to far-IR SEDs of the 43 objects in the 
redshift sample. They are split into 3 luminosity classes - Star- 
bursts (red), LIRGs (green) and ULIRGs (blue). Due to the flux 
density limits of the 70 fim survey, an increase in total infrared lu- 
minosity corresponds to increasing redshift, with all star bursts at 
z < 0.25 and aU ULIRGs at z >0.65. The AGN-dominated SED 
(see section 2.2) is shown in black. The photometry {BRIJK, 3.6, 
4.5, 5.8, 8, 24, 70 and 160 fim) is rcdshift-corrected. 



galaxies, we require models of radiative transfer through 
dusty media, considering solely starburst processes. Ac- 
cordingly, object 165, whose SED we evaluate as AGN- 
dominated (see section 2.2), is not included in our analysis. 
Nevertheless, we do not exclude a strong AGN contribution 
in the remaining galaxies. Simply, at this stage, due to the 
extremely obscured nature of these objects, we cannot reach 
any further conclusions without exploring their multiwave- 
length properties, which is beyond the scope of this paper. 

We use the Chary & Elbaz (2001, hereafter CEOl), 
Dale & Helou (2002, hereafter DH02) and Siebenmorgen & 
Kriigel (2007, hereafter SK07) libraries, selected on account 
of their diverse frameworks and the fact that they span a 
wide range in infrared luminosity (10^ -W^'^Lq), from star- 
forming galaxies to HyLIRGs. In subsection 3.1 we touch 
upon the main characteristics of the models that are rele- 
vant to our study. For a more detailed report, we refer the 
reader to the relevant papers. 

3.1 The dust models 

The CEOl SED templates are based on the models of Silva 
et al. (1998, hereafter S98). S98 employ a stellar population 
synthesis framework, using Simple Stellar populations (SSP) 
of different age and metallicity, including circumstellar dust 
shells/outflows of AGB stars that continuously enrich the 
ISM. Two types of environments are responsible for the in- 
frared continuum - dust emission associated with starform- 
ing molecular clouds and a cirrus component associated with 
the intensity of the stellar field. CEOl generate template 
SEDs, in the UV-submm range, of four prototypical galax- 
ies, Arp 220, NGC 6090, M82, and M51, corresponding to 
four luminosity classes: ULIRGs, LIRGs, starbursts and nor- 
mal galaxies. They partition the four SEDs into mid and far- 
infrared components and interpolate between them to create 
a set of templates with a range of luminosities, from which 
they select the ones which correspond to observed ranges of 
IRAS and Infrared Space Observatory (ISO) colours. 

The DH02 templates are based on the work of Desert, 
Boulanger & Puget (1990, hereafter DBP90). DBP90 in- 
troduce an empirical model combining emission from Very 
Large Grains (VLGs) in thermal equilibrium, stochastically 
heated Very Small Grains (VSGs) and Polycyclic Aromatic 
Hydrocarbon (PAH) features, in radiation fields of vary- 
ing intensity. DH02 produce a set of semi-empirical infrared 
SEDs, the shape of which is determined by the intensity U 
of the interstellar radiation field. They develop a power law 
dependence of the dust mass on the value of U, the outcome 
of which is a collection of dusty environments immersed in 
a range of radiation fields, representing a mixture of active 
and quiescent regions. Their SEDs are further constrained 
by IRAS and ISO observations of local star-forming galax- 
ies with LiR ranging from 10** to 10^^ Lq. In addition they 
vary the emissivity index /3, which relates the emissivity of 
dust grains to the wavelength of emission. 

The SK07 model SEDs are built on the formulation of 
Kriigel & Siebenmorgen (1994), who present a modified def- 
inition of a starburst nucleus: instead of a central radiating 
point source, they employ a configuration where stars are 
distributed throughout the volume, introducing 'Hot Spots' 
- dust regions enveloping OB stars. Correspondingly, in their 
calculations, the dust temperature is varied both in the lo- 
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Figure 3. SEDs for the first 21 objects from the redshift sample fitted with the CEOl (black dashed curve), DH02 (red dotted curve) 
and our empirical templates (blue dot-dashed curve). Only the best-matched (lowest x^) template is shown. The filled circles are Spitzer 
broadband data for the redshift sample, including calibration errors of 10 per cent for the IRAC and MIPS 24 /im photometry, 20 per 
cent for MIPS 70 /im and 30 per cent for MIPS 160 /tm. All photometry has been colour-corrected and rcdshift-corrected so the plots are 
at z~0. 



cality of each star as well as with distance from the nucleus. 
Furthermore, since part of the stellar flux is absorbed within 
the hot spot and re-radiated in the IR, their formulation re- 
places some OB stars by IR-emitting hot spots. This has 
the effect of reducing the overall temperature, whilst keep- 



ing mid-IR emission unchanged, as this is considered to be 
linked to the density of dust within the hot spot. All other 
stars are treated separately and are assumed to contribute 
to the far-IR part, 'cirrus' component of the SED. Their li- 
brary of ~ 7000 templates ranges in 5 free parameters, in 
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Figure 4. SEDs for the remaining 21 objects from tlie redsliift sample fitted witli the CEOl (black dashed curve), DH02 (red dotted 
curve) and our empirical templates (blue dot-dashed curve). Only the best- matched (lowest x^) template is shown. The filled circles are 
Spitzer broadband data for the redshift sample, including calibration errors of 10 per cent for the IRAC and MIPS 24 ^m photometry, 
20 per cent for MIPS 70 ^m and 30 per cent for MIPS 160 fim. All photometry has been colour-corrected and redshift-correeted so the 
plots are at z~0. 



physically acceptable combinations: radius of dust emitting 
region, total luminosity, visual extinction and hot spot dust 
density. The fifth parameter, Lqb / Ltot, represents the per- 
centage of the total infrared luminosity that originates from 
OB stars, compared to the contribution from the general 



stellar population and a result it has a direct effect on the 
near-IR flux and the representation of the near-IR stellar 
bump. 
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Figure 5. SEDs for the first 21 objects from the rcdshift sample fitted with the SK07 templates. The red curves represent the matches 
with the lowest values (see section 3.3), fitted on 8, 24, 70 and 160 /im photometry. The green curve represents the match with the 
lowest value when fitted on 3.6, 4.5, 5.8, 8, 24, 70 and 160 fim photometry. The filled circles are Spitzer broadband data for the 
redshift sample, including calibration errors of 10 per cent for the IRAC and MIPS 24 fim photometry, 20 per cent for MIPS 70 fim and 
30 per cent for MIPS 160 fira. All photometry has been colour-corrected and redshift-corrected so the plots are at z~0. 



3.2 Our empirical templates 

In addition to tiie models and as a consistency check, 
we fit a set of real galaxy templates, which we construct 
by combining spectral and photometric data of local star- 



bursts. The templates are built on a collection of mid-IR (5- 
38 fim) Spitzer InfraRed Spectrograph (IRS) spectra (cour- 
tesy of Brandl et al. 2006) from the central region of 16 
local Starbursts, as well as IRAS broadband data (12, 25, 
60 and 100 fim) from the NASA Extragalactic Database 
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Figure 6. SEDs for the remaining 21 objects from tlie redshift sample fitted witii tlie SK07 templates. Tlie red curves represent tiie 
matclies witii tiie lowest values (see section 3.3), fitted on 8, 24, 70 and 160 fim photometry. The green curve represents the match 
with the lowest value when fitted on 3.6, 4.5, 5.8, 8, 24, 70 and 160 fira photometry. The filled circles are Spitzer broadband data for 
the redshift sample, including calibration errors of 10 per cent for the IRAC and MIPS 24 fira photometry, 20 per cent for MIPS 70 ^m 
and 30 per cent for MIPS 160 fim. All photometry has been colour-corrected and redshift-corrected so the plots are at z ~ 0. 



(NED). The far-IR/submm region (longward of 60 /xm) is 
fitted with a modified black-body emission law, as it is 
considered to sample emission from cool dust in equilib- 
rium heated by the interstellar medium. This is of the form 
B\{T){l — e~^^), with a wavelength dependent optical depth 



Tx = Tioo/jm (lOO/xm/A)'' (e.g. see Klaas et al. 2001) and a 
dust emissivity index (3. We assume a low opacity limit, so 
approximate the term (1 — e"'^-^) by A"''. Typical reported 
values of P range between 1.5-2 (e.g. Dunne et al. 2000, 
Lisenfeld, Isaak & Hills 2000), so we adopt 13=1.5, consistent 
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with studies of the far-IR emissivity of large grains (Desert, 
Boulanger & Puget 1990). The greybody temperature is al- 
lowed to vary in order to minimise on the 60-100 ^m re- 
gion, resulting in a temperature range of 30-50 K. These 16 
empirical SED templates vary with respect to total infrared 
luminosity, position of far-IR peak and mid-IR continuum 
slope. 



3.3 Model Comparison 

The best (lowest x^) SED matches for the redshift sample 
are shown in figures [3Hn] AH photometry has been colour- 
corrected according to the specifications in the IRAC and 
MIPS data handbooks. With the CEOl and DH02 SEDs 
we use all available IRAC photometry and the three MIPS 
bands. In the case of our empirical SEDs, only the 8-160 fim 
region is taken into account, because the templates do not 
extend to short wavelengths. The SK07 models do not treat 
the near-IR emission in detail, so the fits are approached in 
two ways: we fit the 3.6-160 fim region (where IRAC pho- 
tometry is available) and the 8-160 fj,m region separately. 
For the latter all matches with similarly low values are 
kept, in order to increase the accuracy of the results, since 
the library contains ~7000 SEDs and limiting the number 
of photometric points for the x^) dramatically reduces the 
constraints on the chosen templates. Furthermore, near-IR 
photometry is fit with slightly increased errors in order to 
take into account the weakness of the models in that region. 

In the mid-IR (Ares* < 40 fim) all templates deliver the 
emission/absorption lines and continuum slope with the 
same consistency; in some cases the 8 nm colour-corrected 
flux samples a PAH feature (e.g. objects 98, 100, 103, 122, 
125). However, there are large discrepancies in the 10 fim sil- 
icate absorption depth. For the CEOl and DH02 templates 
it seems to vary over a small range; this is not the case for 
the SK07 SEDs, where it can be 50-100 times below the 
continuum (e.g. see objects 58, 70, 76, 82, 110). This is a di- 
rect consequence of SK07 incorporating dust self-absorption 
in their models, an effect assigned to increasing quantities of 
dust at low temperatures (e.g. Mitchell & Robinson 1981), 
resulting in an overall shift of the far-IR peak to longer wave- 
lengths. Our photometry misses the silicate feature so we 
cannot constrain its depth, nevertheless its varied strength 
in the models is consistent with observations as it is com- 
monly found in LIRG & ULIRG spectra (e.g. Hao et al. 
2007, Armus et al. 2007) and its existence is correlated with 
the magnitude of extinction. 

The differences between the templates are more strik- 
ing in the far-IR {Xrest >40/im). The SK07 models are in 
good agreement with the data. In addition, most 3.6-160 /im 
fits result in templates with similar characteristics to the 8- 
160 nm fits, showing that the far-IR is well constrained even 
without the < 6 /im region. On the contrary, DH02 and CEOl 
often underproduce the 160 /im flux, by an amount signifi- 
cantly higher than the 30 per cent photometric uncertainty 
and up to a factor of 10 in extreme cases (e.g. see objects 58, 
70, 84, 93, 94, 110, 134, 137, 138, 152, 172 for which the ef- 
fect is more exaggerated). Moreover, it is worth noting that 
most objects for which the far-IR is completely underesti- 
mated by CEOl and DH02, have a ratio of L7o/i24 > 3.4 
(table 3), suggesting that a steep mid-IR continuum and 



a high far-IR flux are contradictory features in these tem- 
plates. 

Table 2 shows the parameters generated from the best 
matched SK07 templates. The nuclear radius represents the 
region of total dust emission from the entire stellar popu- 
lation, with the OB stars only confined to the inner 350pc. 
The majority of best fit models are described by a radius 
of 3kpc, suggesting emission from extended regions of cool 
dust. In terms of visual extinction, the matches span the 
whole available range of 2 < j4„ < 144. Taking into account 
the variations of A„ for each object with respect to each 
template, it is worth noting that only 5 out of 42 galaxies 
are matched solely with < 10 templates, the remaining 
appearing mostly with Av 35 or 70, implying high obscu- 
ration in these systems. Finally, the value of hot spot dust 
density, which is directly related to the strength of mid-IR 
emission appears to have a wide range in values, but in most 
objects it reaches the maximum of p= 10000 cm"''. 

3.4 Summary 

For the 8-160 fim region, the SK07 models provide the best 
matches to our sources' infrared SEDs. The fact that the 
CEOl and DH02 templates in most cases severely underes- 
timate the far-IR, leads us to two conclusions. Firstly, the 
notion of a single radiating source is an over-simpliflcation 
of a real system, which is more likely to be a complex and 
clumpy 3-D conflguration with a varying temperature gra- 
dient. Secondly, we propose that these high-z objects have 
higher obscuration and additional far-IR emission from a 
cold dust component that is not commonly seen in local 
LIRGs/ULIRGs, on which the CEOl and DH02 templates 
are based. Such excess far-IR flux would be a direct conse- 
quence of a large emitting dust region, where dust in the 
immediate vicinity of young stars (and/or an AGN), would 
be responsible for greatly attenuating the radiation held and 
hence shielding the outer dust layers from direct UV illu- 
mination, resulting in lower equilibrium temperatures. This 
scenario is well-expressed by the SK07 models; by includ- 
ing the treatment of hot spots around OB stars, the aver- 
age dust temperature is lowered, shifting the far-IR peak 
to longer wavelengths, without a reduction in the mid-IR 
flux, which, according to their framework, is only linked to 
the dust density in the hot spot. In contrast, the CEOl and 
DH02 templates, which have been found to work well with 
local sources, fail in this case, because their formulation does 
not allow for the possibility of excess far-IR flux independent 
from the mid-IR. These observations are further supported 
by our empirical SED fits, which were included for the pur- 
pose of comparing with real local galaxy data (see section 
3.2). The majority show considerable discrepancies, failing 
to reproduce the far-IR SED shape (e.g. objects 77, 93, 110, 
122), indicating that the greybody temperatures of 30-50 K 
are too high and that dust potentially colder than 30 K is 
present in most systems. 

In view of these results, our subsequent calculations 
consider only the SK07 fits to the 8-160 /im spectral re- 
gion. The fact that the near-IR part of the SED is not well- 
represented, is not important since we are only concerned 
with pure dust emission in the 8-1000 fim range. For each 
object, we take into account all SK07 matches (figures [5] 
and [6] and table 2), to establish a valid error margin in cases 
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Table 2. Parameters corresponding to the SK07 templates, fitted with photometry at 8, 24, 70 
and 160 jum, that produced the lowest matches (used in subsequent calculations). The columns 
are as follows: Object ID, radius of the central nuclear region (kpc), ratio of the luminosity from 
OB stars to the total luminosity (LoB/Ltot), total visual extinction measured from the surface 

to the centre of the starburst (Ay) and density within the hot spots p (cm~^). Where a range is 
quoted, e.g. p=100— 10000, each of the matched models came up with a value of p between those 
limits. Where there is a coma, e.g. 100, 10000, each of the matched models had either p of 100 or 
10000, but not in between. 



ID 


IN LiCiedl I dUl Uo 


^OB / i^tot 


A 


P 




l^Kpc; 


^/o; 




(cm ) 


41 


0.35,3 


40,60 


35.9,72 


100-5000 


51 


3 


40,90 


35.9,67.3 


10000 


55 


0.35,3 


40,60 


35.9,72 


100-10000 


58 


0.35,1,3 


40,60 


119,144 


5000,10000 


67 


1,3 


60 


35.4,35.9 


10000 


70 


0.35,1 


40,60 


119,144 


5000,10000 


72 


1,3 


40 


2.2-35.9 


1000-10000 


76 


0.35,1 


40,60,90 


35.9-119 


100,10000 


77 


1,3 


40,60 


35.9-72 


2500, 10000 


78 


3 


40,60 


35.9,72 


2500,10000 


82 


1,3 


60,90 


18-72 


100-10000 


84 


1,3 


60 


70.7,72 


10000 


85 


1,3 


40 


2.2,35.9 


1000,10000 


88 


0.35,3 


60,90 


72,120 


10000 


93 


0.35,1 


40,60 


35.4-72 


10000 


94 


1,3 


40 


35.4,35.9 


10000 


98 


1,3 


40,60 


4.5-35.4 


2500,10000 


100 


3 


60,90 


2.2-6.7 


100 


103 


1,3 


40,60,90 


2.2-70.7 


100-10000 


110 


0.35,1 


40,60,90 


35.4-144 


1000,10000 


114 


3 


40 


17.9 


2500 


115 


3 


40,60 


2.2-72 


1000,10000 


120 


1,3 


40 


2.2-9 


1000,2500 


121 


1,3 


40,60 


9-35.9 


5000 


122 


0.35,3 


40,60 


35.9-120 


5000,10000 


124 


1,3 


40 


17.9,35.9 


7500,10000 


125 


3 


60 


2.2,4.5 


100 


126 


0.35,1 


40,60,90 


4.5-17.9 


100 


134 


1,3 


40,60 


70.7,72 


10000 


137 


1,3 


40 


17.9-72 


1000,10000 


138 


1,3 


40 


35.9-72 


1000,10000 


140 


1,3 


40,60 


9-70.7 


1000-10000 


141 


0.35,1 


40,60 


2.2-35.9 


100-10000 


144 


3 


40,60 


2.2 


100 


146 


1,3 


40 


4.5-35.9 


1000,10000 


147 


3 


40 


17.9,35.9 


2500 


149 


1,3 


40 


6.7-35.9 


1000-10000 


150 


1,3 


60 


17.9,35.9 


7500,10000 


152 


1,3 


40,60 


70.7,72 


10000 


155 


0.35,1 


40,60,90 


2.2-6.7 


100 


158 


0.35,3 


40,60 


35.9,72 


2500-10000 


172 


0.35 


40,60,90 


144 


2500 



where the 160 fim detection could either reside on the peak 
or form part of the rising continuum (e.g. see objects 115, 
147, 152). 

4 SED CHARACTERISATION 

4.1 The Infrared Energy Budget 

With the aid of the best-matched SK07 templates, we per- 
form K-corrections and estimate the monochromatic rest 
frame luminosities at 24, 70 and 160 /xm. In addition, by 



integrating the SEDs between 8-1000 nm, wc arrive at an 
array of values for the total infrared luminosity {Lm) for 
each object, the average of which is then used as a final es- 
timate. The sources in the rcdshift sample arc categorised 
into luminosity classes, following the definition in Sanders & 
Mirabel (1996), revealing that it consists of 12 per cent Star- 
bursts, 62 per cent LIRGs and 26 per cent ULIRGs. Despite 
the uncertainties on Lir being relatively large and stem- 
ming from the lack of photometry longward of 160 ^m, only 
two objects move to a different luminosity class on account 
of this (table 3) . 
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Figure 7. The K-corrected, rest frame luminosity ratios 1/70/^24, iieo/^TO as a function of the total infrared luminosity. All units are 
Lq and the Icr uncertainties are shown. 
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Figure 8. The position of the rest frame SED peak (in fim) as a function of the K-corrected, rest frame L7o/i24 ratio (left panel) and 
Lifio/Lro ratio (right panel). L24, L70 and Lieo are in units of Lq and the Icr uncertainties are shown. The colours correspond to the 
luminosity class: black for Starbursts, red for LIRGs and green for ULIRGS. 



4.2 The mid-IR continuum slope and SED peak 

Conventional methods of characterising infrared galaxies in- 
clude associating flux ratios between various bands to phys- 
ical processes responsible for emission in the infrared (e.g. 
Xu et al. 2001; Egami et al. 2004; Lutz et al. 2005; Verma et 
al. 2005). In particular, the /25//60 parameter has been ap- 
plied extensively in IRAS studies of such systems, to quan- 
tify relative AGN/starburst contributions with an indica- 
tive value of /25//60 > 0.17 for AGN-dominated sources (e.g. 
de Grijp et al. 1985, Miley et al. 1985, Soifer et al. 1989). 
This 'warm/cold' terminology is based on the implication 
that dust heated by an AGN would reach higher tempera- 
tures and hence lower the continuum slope between 25 and 
60 pm (or 24 and 70 fj,m in Spitzer^s case) (e.g. Farrah et al. 
2005; Verma et al. 2005, Frayer et al. 2006). Calculating the 
25/60 /im continuum slope using the model SEDs, reveals 
that for ~25 per cent of our sample the classical IRAS warm 
source criterion of /25//60 > 0.17 apphes. However, the ma- 
jority of these objects peak at A > 100 nm, indicating that 
the flatter behaviour between 25 and 60 /im is due to a shift 
in infrared emission with elevated A > 90 /im flux and re- 
duced 60 /im flux. Furthermore, from figure [T] it is clear that 
the continuum slope throughout the SED is completely in- 
dependent from the total infrared luminosity. Consequently, 
if it is true that the contribution of AGN-dominated ob- 
jects increases as a function of IR luminosity (e.g. Brand 
et al. 2006), then the slope of the mid-IR continuum and 



hence the L70/L24 ratio, fail to quantify it. In addition, be- 
cause of heavy obscuration, emission from cool dust seems 
to be altering the mid-IR continuum slope solely due to in- 
creased far-IR emission. Various authors have supported this 
idea; Chakrabarti et al. (2007) have modelled infrared SEDs 
with varied starburst/ AGN contributions, suggesting that 
the warm/cold classification of LIRGs/ULIRGs would be 
independent of the energy source, if sufficiently obscured so 
that all the far-IR emission comes from large grains in equi- 
librium. Similar observations (e.g. Lutz et al. 1996, Klaas et 
al. 1997) have put forward the notion that elevated optical 
extinction can be responsible for a steep rise of the mid- 
IR continuum and hence can displace the ratio of the 'cold 
-warm' components, resulting in an energy output mainly 
carried by the cold component. 

To further investigate the properties of this sample, the 
wavelength at which the SEDs turn over into the Rayleigh- 
Jean regime is estimated (table 3, column 8). We find that 
a significant fraction of objects (40-60 per cent within Icr) 
peak at wavelengths longward of 90 fim independent of lumi- 
nosity class. Despite the large uncertainties, there is clearly 
a strong correlation with respect to I/i6o/i70, not surpris- 
ing, since 70 and 160 fim sample either side of the peak and 
their ratio is directly related to the far-IR turnover (figure 
[SI right panel). The same does not apply for the L70/I/24 
ratio, however, which shows no such dependence (figure |8l 
left panel). 
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Table 3. The infrared budget and SED characteristics of the redshift sample. Column 1 classifies 
the objects according to their luminosity. Columns 2 and 3 are the Object ID and redshift, columns 
4 and 5 are the K-corrected, rest-frame L70/L24 and Lieo/Lro ratios. Columns 6 and 7 quote 
the average value of the total infrared luminosity Ljji (in units of xIO^'^Lq) and its logarithm. 
The position of the far-IR turnover is displayed in column 8. All uncertainties are Icr. 



Type 


ID 


L7o/L2i 




LiR X 1O1O(L0) 


logLiR 


peak (rest, fim) 




51 


1.86 


1.52 


4.32± 0.23 


10, 


.64 


140± 14 




155 


3.00 


0.17 


2.89± 0.2 


10, 


.46 


55± 2 


SBs 


125 


4.10 


0.55 


6.35± 0.72 


10, 


.81 


79± 5 




144 


5.56 


0.61 


4.10± 0.49 


10, 


.61 


86± 2 


SB /LIRG 


103 


3.45 


0.95 


8.79± 1.31 


10, 


.94 


98± 28 




141 


0.96 


0.10 


36.54± 1.66 


1 1 
11 


. 00 


41± 3 




121 


1.59 


1.04 


72.07± 8.26 


1 1 
11 


.OD 


111± 17 




150 


1.62 


1.29 


58.88± 9.01 


Y\ 


77 


117± 22 




98 


2.27 


0.75 


11.97± 1.0 




.uo 


95± 7 




100 


2.51 


0.42 


22.68± 1.48 




.ou 


64± 5 




78 


2.62 


2.01 


34.27± 1.03 




.00 


141± 13 




120 


2.80 


0.49 


31.15± 2.84 


Y\ 


AQ 


78± 5 




152 


2.84 


2.09 


37.41± 8.58 


11, 


.57 


140± 18 




115 


2.88 


1.22 


26.50± 7.93 


11, 


.42 


114± 23 




85 


2.88 


0.57 


23.73± 2.32 


^\ 


.37 


84± 8 




76 


3.00 


0.11 


35.41± 3.65 


11, 


.55 


50± 3 




72 


3.02 


0.47 


13.81± 0.87 


11, 


.14 


78± 5 




147 


3.03 


1.26 


74.57± 16.67 


11, 


.87 


113± 22 


LIRGs 


122 


3.24 


1.15 


12. lilt 2.31 


11, 


.08 


109± 13 




77 


3.33 


2.08 


81.11± 16.14 


11, 


.91 


135± 15 




146 


3.34 


0.63 


21.24± 0.98 


11, 


.33 


90± 7 




140 


3.35 


1.29 


31.18± 4.74 


11, 


.49 


119± 16 




149 


3.48 


0.95 


47.50± 7.27 


11, 


.68 


107± 9 




55 


3.59 


0.42 


86.61± 3.77 


11, 


.94 


79± 10 




126 


3.60 


0.23 


30.87± 1.90 


11, 


.49 


58± 2 




41 


3.66 


0.81 


77.76± 8.38 


11, 


.89 


llOit 26 




137 


4.55 


1.41 


63.65± 4.0 


11 


.80 


121± 8 




138 


6.27 


0.77 


15.82± 0.59 


11, 


.20 


99± 6 




110 


7.58 


0.25 


33.68± 3.37 


11, 


.53 


67± 3 




172 


13.69 


0.50 


37.65± 0.12 


11, 


.58 


85± 


(U)LIRG 


134 


3.79 


1.48 


83.94± 21.39 


11, 


.92 


119± 15 




67 


2.58 


0.56 


224.71± 2.74 


12, 


.35 


86± 9 




114 


2.71 


1.27 


135. 64± 1.54 


12, 


.13 


119± 9 




82 


2.86 


0.30 


283. 69± 39.97 


12, 


.45 


57± 3 




124 


2.89 


0.38 


146. 39± 8.87 


12, 


.17 


75± 5 


ULIRGs 


93 


3.04 


0.40 


493. 36± 55.53 


12, 


.69 


76± 7 




84 


3.43 


1.27 


202. 19± 8.25 


12, 


.31 


118± 13 




88 


3.77 


0.23 


242. 27± 15.29 


12, 


.38 


65± 2 




94 


4.03 


0.49 


273. 44± 2.96 


12, 


.44 


82± 5 




158 


4.71 


0.43 


506. 70± 111.65 


12, 


.70 


77± 20 




70 


9.40 


0.69 


201. 46± 27.19 


12, 


.30 


93± 15 




58 


10.48 


1.06 


121. 73± 16.74 


12, 


.08 


108± 18 



Figures [7] and [8] display some important properties of 
the sample, consistent with conclusions from section 3. Gen- 
erally, a simple black body spectrum displays a temperature 
dependence through the position of the emission peak, the 
continuum slope and the total luminosity. Consequently, if 
the far-IR part of an SED were well-represented by a [mod- 
ified] black body, for a sample of galaxies of varying dust 
temperatures this would manifest itself as: a) a positive cor- 
relation between mid-IR continuum slope and total infrared 
luminosity (as they both increase with temperature) and b) 
a negative correlation between mid-IR continuum slope and 
position of SED peak (as the former increases and the lat- 



ter decreases with temperature). However, if an extra cold 
dust component is present, as we believe the case to be for 
the 70 /im sample, a single temperature black body is not a 
sufficiently accurate representation of the far-IR and these 
correlations would not apply. This is what is observed in 
figures [7] and |8l 



5 EQUATIONS FOR ESTIMATING THE 
TOTAL INFRARED LUMINOSITY 

Determining the infrared energy budget is crucial for esti- 
mating the contribution to the infrared and submillimetre 
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Figure 9. The relation between Lju calculated using equation 
1 (y-axis) and Lm calculated by integrating under the SEDs (x- 
axis). The solid line has a gradient of 1. The r.m.s. scatter off the 
line is 0.06 dex. 

backgrounds and plays an important role in establishing a 
value for the infrared element of the cosmic star formation 
history. If the available photometric data is not sufficient 
and hence unable to explicitly define these contributions, a 
first order approximation of the infrared energy budget can 
be attained through equations relating the total infrared lu- 
minosity to the monochromatic luminosity in one or more 
bands. In this section, we develop equations that can be 
used to estimate Ljr, using 1 or more MIPS bands. How- 
ever, note that all have been derived using the typical form 
of starburst-type SEDs and have not been tested on power- 
law SEDs of AGN-dominated sources. 



5.1 Estimating Lir with all 3 MIPS bands 

As the SED of an infrared-luminous galaxy has a distinct 
shape, its bolometric output can be approximated by sum- 
ming contributions from four main sections: < 24 ^m, 24- 
70 ^m, 70-160 /im and > 160 /im. In a linear plot of fv ver- 
sus f, flux shortward of 24 nm and longward of 160 /im is 
not as significant as the fiux within that range. Therefore 
the < 24 /im and > 160 /im sections can be approximated by 
triangles of heights L24 and Lieo and the sections within 24- 
160 /im with trapezia of heights of L24, L70 and L70, I/ieo- 
Summing these contributions and empirically accounting for 
the systematic offset in this approximation, we derive an 
equation to calculate the total infrared luminosity {Lir) in 
the range 8-1000 ^m: 

Lir = 4.63 X 10~'^ x (8.3L24 + 2.7L70 + Lieo) x (1 + z).{l) 

The coefficients represent the width (Ai^) of the bases of 
the sections in units of Hz, Lir is in units of Lq and L24, 
L70 and L160 are the monochromatic luminosity densities 
(in units of W/Hz) directly derived from the observed fiux 
densities, without applying K-corrections. The factor (1-l-z) 
adjusts the bases Av of the SED sections, as their width 
increases with redshift. 

The validity of this method is demonstrated by compar- 
ing the value of Lir derived by using the equation (y-axis) 
and the average value of Lir derived by integrating under 
the best-matched model SEDs (x-axis) (figure There is 
very good agreement, with an r.m.s. deviation of 0.06 dex 
from the true Lir. We also test it on a different galaxy 




log L|R (Lq) 



Figure 10. The relation between Lir calculated using equation 
1 (y-axis) and Lir calculated by integrating under the SEDs (x- 
axis) for 21 SINGS galaxies. The solid line has a gradient of 1. 
The r.m.s. scatter off the line is 0.05 dex. 



sample, by calculating the total infrared luminosity of 21 
local infrared star-forming/starburst galaxies using MIPS 
photometry from the Spitzer Infrared Nearby Galaxy Sur- 
vey (SINGS) (Dale at al. 2005, hereafter DOS) and estimated 
distances quoted in Sanders et al. (2003). We fit the DH02 
models on the 24, 70 and 160 fim photometry (like in DOS) 
and integrate under the 8-1000 /xm region of the best-fit 
curve, comparing the result with Lir from equation 1. The 
agreement is very good, with an r.m.s. scatter of O.OS (figure 
nop , confirming that this equation is not specific to high-z 
LIRGs/ULIRGs, but also applicable to normal infrared star- 
forming galaxies in the local universe. 



5.2 Estimating Lir with 1 MIPS band 

Naturally, a more luminous galaxy will also display higher 
luminosity at the various infrared photometric bands (e.g. 
Elbaz et al. 2002, Takeuchi et al. 200S). This would mani- 
fest itself as a strong positive correlation between Lir and 
K-corrected rest-frame monochromatic luminosity at 24, 70 
and 160 ^im (figure [TT|) . As a result, it could be used instead 
of equation 1, for estimating Lir, in cases where only one 
photometric data point from MIPS is available. For each 
plot in figure [TT] we compute the least squares best fit line 
through the data, equations 2-4: 

log Lir = 1.9S + 0.89 log L24, (2) 

log Lir = 1.16 + 0.92 logLro, (3) 

logLiR = 1.49 + 0.9 log L160, (4) 

All luminosities are in units of Lq and L24, I/70 and Lieo are 
in the rest frame. The r.m.s. offsets from the best line fits 
are 0.18, 0.11, 0.23 dex for equations 2, 3 and 4 respectively. 
This intrinsic dispersion most likely stems from the varia- 
tion in the properties of each galaxy, the high photometric 
uncertainties as well as the range of models used to calcu- 
late each luminosity value; nevertheless, they are convenient 
in cases where the luminosity of a galaxy is required only 
within an order of magnitude or class. 
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Figure 11. The monochromatic rest frame luminosity at 24, 70 
and 160 fim as a function of the total infrared luminosity for the 
redshift sample. The Icr uncertainties are shown. The r.m.s. scat- 
ter is 0.18, 0.11, and 0.23 dex respectively. 



6 SUMMARY AND CONCLUSION 

We have examined a population of 70 ^im selected sources, 
in the region of the Extended Groth Strip, focusing this 
study on the infrared properties of 43 objects with available 
spectroscopic redshifts (0.1 < z < 1.2). In the last section, we 
derived a set of equations in order to estimate the total IR 
luminosity with at least one MIPS band, applicable to ob- 
jects with starburst-type SEDs. We would like to emphasise 
that our selection criteria did not photometrically bias the 
redshift sample with respect to the full sample or predis- 
pose our study to a particular type of galaxy. As this is a 
flux limited survey (see section 2), the higher redshift objects 
are also more luminous in the infrared and display lower op- 
tical flux. This implies that the sources with no available 
spectroscopic redshift are also in the LIRG/ULIRG regime. 
In addition, the MIPS colours (figure [l| and SEDs (figure 
[5| unanimously show that all objects are photometrically 
similar in the mid and far-IR and their properties diverge 
only in the optical and near-IR, not relevant for the work 
described here. 

By fitting 4 Ubraries of SED templates on our infrared 
photometry, we have examined the properties of the sources 
and estimated the infrared energy budget, revealing 12 per 
cent starbursts, 62 per cent LIRGs and 26 per cent ULIRGs. 
Evaluation of the fits showed that the SK07 models perform 



better at reproducing the far-IR SED; their formulation 
adopts a two-component configuration of dust heated in the 
immediate locality of massive young stars and cirrus-heated 
dust, essentially decoupling mid and far infrared emission. 
As a result, since the total infrared luminosity is not directly 
related to the position of the SED peak or the steepness of 
the mid-IR continuum slope, with the SK07 templates it is 
possible to represent objects with strong mid-IR emission 
and an additional cold component. CEOl and DH02 em- 
ploy a scenario with a single radiating central source where 
temperature varies as a function of distance from the centre. 
Although this has worked well with local infrared galaxies, it 
completely underestimates our sources' 160 /im fiux, with a 
discrepancy of at least a factor of 2 and up to 10 for extreme 
cases, since it does not allow for types of systems with ele- 
vated far-IR emission. Our empirical local galaxy templates, 
also mostly underestimate the far-IR region, indicative of 
the fact that an average dust temperature of 30-50 K is too 
high. 



Considering all above points, we propose that this deep 
70 /im survey has probed high-z LIRGs and ULIRGs, which 
are unlike local counterparts: heavy obscuration and large 
amounts of cold dust potentially lower than 30 K, appear 
as a far-infrared excess component in the SEDs, causing 
a significant fraction to peak at A > 90 /im and steepen- 
ing the continuum slope. The existence of such objects has 
also been put forward by other authors; e.g. Marcillac et al. 
(2006) have suggested that some high-z infrared-luminous 
objects display far-IR properties divergent from those of lo- 
cal galaxies because of an additional cold dust component, 
even if the strength of mid-IR emission is the same. This 
has further implications: in such systems, the continuum 
slope cannot be representative of various galaxy character- 
istics such as dust temperature, total infrared luminosity 
or the nature of the central energy source. Figures [7] and [8] 
demonstrate this: the mid-IR continuum slope (quantified by 
the Z/70 / Z/24 ratio) satisfies the IRAS warm source criterion 
of /25//60 > 0.17 for sources with reduced 60 /im emission 
and increased A > 90 fim fiux, as opposed to elevated 25 /im 
emission and it is completely dissociated from the position 
of the SED turnover and total infrared luminosity. As this 
is inconsistent with the representation of far-IR emission as 
a single temperature black body, it is in line with our ear- 
lier suggestion of the presence of an additional cold emissive 
component. 



It is worth noting that this 70 /im sample is one of the 
first far-IR-selected samples reaching such low fiux-density 
limits and additional similar surveys will be of great impor- 
tance in highlighting the differences (if any) between local 
and high redshift infrared populations. The next infrared ob- 
servatory, Herschel, will greatly benefit this work, extending 
to the sub-mm part of the spectrum and, hence, enabling a 
complete census of high-z infrared galaxy SEDs. Finally, we 
would like to mention that in order to confirm our conclu- 
sions we have planned a follow-up of the work described in 
this paper with submillimetre observations using the Sub- 
millimetre Common-User Bolometer Array (SCUBA 2). 
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